Analytical ground state for the Jaynes-Cummings model with the ultrastrong coupling 
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We present a generalized variational method to analytically obtain the ground-state properties 
of the unsolvable Jaynes-Cummings model with the ultrastrong coupling. An explicit expression 
for the ground-state energy, which agrees well with the numerical simulation in a wide range of the 
experimental parameters, is given. In particular, the introduced method can successfully solve this 
Jaynes-Cummings model with the positive detuning (the atomic resonant level is larger than the 
photon frequency), which can not be treated in the adiabatical approximation and the generalized 
rotating-wave approximation. Finally, we also demonstrate analytically how to control the mean 
photon number by means of the current experimental parameters including the photon frequency, 
the coupling strength, and especially the atomic resonant level. 
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The Jaynes-Cummings model, which describes the im- 
portant interaction between the atom and the photon 
of a quantized electromagnetic field, is a fundamental 
model in quantum optics and condensed-matter physics 
as well as in quantum information science. In the optical 
cavity quantum electrodynamics, the atom-photon cou- 
pling strength is far smaller than the photon frequency. 
As a result, the system dynamics can be well governed 
by the Jaynes-Cummings model with the rotating-wave 
approximation (RWA) [1|. In the case of the RWA, its 
energy spectrum and wavefunctions can be solved exactly 
[2| . With the rapid development of fabricated technique 
in solid-state systems, the Jaynes-Cummings model can 
be realized in semiconductingdots [3H6| and supercon- 
ducting Josephson junctions [7Hl2|. More importantly, 
recent experiment has reported the existence of the ul- 
trastrong coupling with the ratio 0.12 between the cou- 
pling strength and the microwave photon frequency [l3| ■ 
Moreover, this ratio maybe approach unit due to the cur- 
rent efforts HI Ell . 

However, in this ultrastrong coupling regime the well- 
known RWA breaks down and the whole Hamiltonian is 
written as 

H x = uj{a)a + -) + -n<j z + g{a+ + ( r_)(a t + a), (1) 

where a^ and a are creation and annihilation operators 
for photon with frequency u, er± are the raising and low- 
ering operators of the two-level atom in the basis of a z , 
17 is the atomic resonant frequency and g is the atom- 
photon coupling strength. Due to the existence of the 
counter- rotating terms (c+a + and cr_a), Hamiltonian ([1]) 
is very difficult to be solved analytically except for O = 0. 
Although the energy spectrum of Hamiltonian ((TJ) has 
been obtained perfectly by means of the numerical simu- 



lation |16l - tT8| , the analytical solutions are very necessary 
for extracting the fundament al phys ics as well as in pro- 
cessing quantum information |19l421| . In the negative de- 
tuning 17 < u, the adiabatic approximation method that 
the second term of Hamiltonian (QJ) is treated as a small 
perturbation has been considered [22|. Recently, a gen- 
eralized rotating-wave approximation (GRWA) has also 
been proposed to solve Hamiltonian (QJ) in the displaced 
oscillator basis states [23|. This method can derive an 
analytical energy of Hamiltonian ([1} in the ultrastrong 
coupling successfully. However, it, like the adiabatic ap- 
proximation method, is invalid in the case of the positive 
detuning (17 > uj) [2J. Moreover, the GRWA also leads 
to an unphysical relation that the mean photon number 
is independent of the atomic resonant frequency. Thus, 
it is very necessary to put forward a new method to re- 
consider Hamiltonian ([TJ. 

In this brief paper we present a generalized variational 
method (GVM) to analytically obtain the ground-state 
properties of Hamiltonian (UJ) in the ultrastrong coupling 
regime [13j . An explicit expression for the ground-state 
energy, which agrees well with the numerical simulation 
in a wide range of the experimental parameters, is given. 
More importantly, our method is valid for all regions of 
the atomic resonant level including the negative detuning 
(17 < w), the resonant case (17 = uj), and especially the 
positive detuning (17 > lu). Finally, we also demonstrate 
analytically that the mean photon number is strongly 
dependent on all parameters of Hamiltonian (QJ) includ- 
ing the photon frequency, the coupling strength, and the 
atomic resonant level. For a weak atomic resonant fre- 
quency, the mean photon number depends linearly on it. 

We first perform a rotation around the y-axis to rewrite 
Hamiltonian (UJ) as H z = uj{a^ a + ^) + ^a x — g<r z {a) +a). 
Under a unitary transformation U — exp[Acr z (a/ — a)], the 



above Hamiltonian H z becomes 
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H u = \ 2 LJ 

with H a — Ld(a^a + | 



2Xg 



-(g + Aw) (at + 



a <t z 



(2) 
and 



H Q = ^Cl{a + exp[2A(a t - a)] + ct_ exp[-2A(at - a)]}, 
where A is a dimensionless parameter. In the Jaynes- 
Cummings model with the KWA, the dimensionless pa- 
rameter A can be determined by setting H g = and an 
effective Hamiltonian can be obtained by series expansion 
with respect to g/ui |7|. However, in the ultrastrong cou- 
pling, this derivation is not valid. Moreover, the terms 
H g and Hn play the important role in the energy spec- 
trum. As will be shown, the dimensionless parameter A in 
the ultrastrong coupling can be determined analytically. 
We now choose the basis states |d=, JV) = 
-4= (c| ± c g ) \0) a \N), where c\ and c g are the cre- 
ation operators of the excited- and ground- states, 
|0) a is the vacuum state of atom and |JV) is the Fock 
state of the photon, to calculate the ener gy spectrum 
by means of the perturbation method [25|]. This 
method requires the diagonal and non-diagonal ma- 
trix elements of Hamiltonian @ in the basis states. 
Since -{g + Aw) (JV, ±| (a f + a)a z \N,±) = and 
u{N,±\(ata + i) \N,±) = uj(N + |), the diagonal part 
is given by 

H a = \ 2 uj + 2\g + uj(N+l) + (Hn) (3) 



with (Hn) 
where F(X) 



(N,±\Hn\N,±) = tF(\)L%(4\ 2 ), 
-±£lexp(-2A 2 ) and L{(x) is the 
associated Laguerre polynomials. The non-diagonal 
parts must be divided into two cases including 
(N,±\H g \±,M) and (N,±\H Q \±,M) for even (JV - 
M), and (N,±\H g \T, M) and (JV, ±| H n |=F, Af) for odd 
(N-M) since (N,±\H a \±,M) = (N 7 ±\H a \t,M) = 0. 
When (JV - M) is even, we have (N,±\H g \±,M) = 
and (N,±\Hn\±,M) = P N>M , where P N<M = 
TF(\) (2\f- M y/M\/N\L^- M (4A 2 ). When (JV - M) 
is odd, the matrix elements become (JV, ±| H g |=p, M) = 
-VNlj( 9 /l; + \)6 n ,m+i and (JV,±| Hn |=F, M) = 
—Pn,m- From the above calculation we assume JV > M 
in terms of the Hermitian of Hamiltonian @. 

By means of the matrix form of Hamiltonian ^ with 
the derived diagonal and non-diagonal terms, the en- 
ergy spectrum and wavefunction can, in principle, be ob- 
tained. However, the procedure is very complicated. For- 
tunately, in the ultrastrong coupling regime, the ground- 
state properties can be calculated by the perturbation 
expansion of the non-diagonal terms. In this method, 
the diagonal part is regarded as the unperturbed Hamil- 
tonian Hq and the remaining part of Hamiltonian, H r = 
H u — H is treated as a perturbation |25| . 

Formally, the ground-state energy can be written as 
261 
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FIG. 1: (Color online) The ground-state energy Eo as a func- 
tion of the atomic resonant frequency Q, for the different cou- 
pling strength g — 0.2w (a) and 0.6w (b) by means of the 
numerical simulation (NS), the GVM and the GRWA. 



The first-order correction to the unperturbed ground- 
state energy is given by Eq = since H r has no diag- 
onal matrix element. While the second-order correction 



is evaluated as E. 



(2) 




(g + u\)+2\F(\)Y/[w 



E 



(0)1 



2F(X) (1 - 2A 2 )]-£~ =2 F 2 (A) (2A) 2iV /&.&%- „ 
where e = ±JV is odd (even). On the other hand, accord- 
ing to the spirit of the variational method (25[, the di- 
mensionless parameter A can be derived from minimizing 
the unperturbed ground-state energy Eq by 



A[cj + £lcxp(-2A 2 )] + <? = 0. 



(5) 



By using the same procedure, the ground-state wave- 
function is written by [261 ] 
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where the unperturbed ground-state wavefunction is 



given by 
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,0), and the first-order correc- 



tion to the unperturbed ground-state wavefunction 



is 



by 
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Equations ((U and ([6]), together with Eq. ((5|), are the 
main results in this brief paper. Interestingly, the di- 
mensionless parameter A depends on all parameters in- 
cluding the photon frequency u>, the coupling strength 
g, and especially the atomic resonant frequency 51. In 
general, Eq. ([5]) can not be solved and the ground-state 
properties can not be obtained analytically due to the 



Here E^ is the unperturbed ground-state energy. Ac- ^™ce the complex sums in the expressions E and 
cording to Eq. §3§ , this unperturbed ground-state energy 
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However, in current experimental setups of the 



is given by Eq — \uj + \ 2 uj + 2A<? + F{\) for JV = 0. ultrastrong coupling (for example, g ~ O.lw in Ref. 1 
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FIG. 2: (Color online) The ground-state energy Eo as a func- 
tion of the coupling strength g for the different atomic reso- 
nant frequency Q. = l.Ow (a), and 1.5w (b) by means of the 
numerical simulation (NS), the GVM, and the GRWA. 



FIG. 3: (Color online) The ground-state energy Eo as a func- 
tion of the atomic resonant level Q. with the coupling strength 
g — I.Olj by means of the numerical simulation (NS), the 
GVM, and the GRWA. 



Eq. ([3]) becomes very simple, and moreover, the complex 
sums are also simplified successfully. As a consequence, 
some interesting results, which govern the fundamental 
properties of Hamiltonian ((T|), can be achieved. 

By means of Eq. ([5]) , the dimensionless parameter A is 
derived approximately by A ~ — g/{uj + ft exp[— 2g 2 /(uj + 
ft) 2 ]}. However, in current experimental setups of the 
ultrastrong coupling, exp[— 2g 2 /(ui + ft) 2 ] ~ 1, and thus, 
the dimensionless parameter A becomes 



(7) 
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Due to the 



existence of the high order terms with respect to A, the 

(2) 

complex sum in Eq can be omitted. Moreover, the 
first term of Eq is far smaller than Eq . It means 
that the second-order perturbed energy almost has no 
effect to the ground-state energy. Therefore, an explicit 
ground-state energy for Hamiltonian ([T]) is given finally 
by 



Eq 



u> g 2 (u + 2Q) 
2 ~ 



n 



{u + ny 



exp[-2(- 



ft' 



(8) 



For a weak atomic resonant frequency (ft <C to), Eq. © 
become E Q 3RWA = l/2w - g 2 /uj - l/2ftexp[-2 (g/u) 2 ], 
which has been obtained perfectly by means of the 
GRWA H|. 

Figure 1 shows the ground-state energy E as a func- 
tion of the atomic resonant frequency ft for the differ- 
ent coupling strength g = 0.2uj and 0.6w by means of 
the numerical simulation, the GVM (Eq) and the GRWA 
(Eq RWA ). For a weak atomic resonant frequency (ft < 



to), both results derived from the GVM and the GRWA 
agree with the direct numerical simulation. In the pos- 
itive detuning (ft > oj), the GRWA breaks down and 
the error is increased linearly. However, the result de- 
rived from the GVM can also agree well with that of the 
numerical simulation. Figure 2 shows the ground-state 
energy Eq as a function of the coupling strength g for 
different atomic resonant frequency ft = l.Ow and 1.5w 
by means of the numerical simulation, the GVM (Eq) 
and the GRWA (E JRWA ). This figure exhibits that the 
analytical expression in Eq. ([HJ can work reasonably for 
g < 0.8cj. If g ~ w, Eq. (O can not be solved explic- 

(2) 

itly and the complex sums in the expressions Eq and 
^g ) can also not be simplified. However, the results 

in Eqs. ([¥]) and ^ with the complex sums remain valid, 
as shown in Fig. 3. 

Having obtained the ground-state energy, we now 
discuss the experimentally-measurable mean photon 
number (a 1 * a), which can be obtained by (a 1 * a) = 
(* | Ua^aW |* > = (*ol o)a + A 2 - \a z ( t + a) |* ). 
The expression of the mean photon number is very com- 
plex. However, here we give an approximate solution, 



(a^a) 



9 



ftexp(-2 5 2 /u; 2 )] 2 ' 



(9) 



Equation ([9]) demonstrates clearly that the mean photon 
number depends strongly on all parameters of Hamilto- 
nian (UJ. For a weak atomic resonant frequency (fi « w), 
the mean photon number becomes (a^a) = g 2 /ui 2 — 
2g 2 flexp(—2g 2 /uj 2 )/uj 3 , which is quite different from the 
previous conclusion (fl'fljg™,, = g 2 /^ 2 derived from 
the GRWA [UJ. 

Figures 4 shows the mean photon number (0,^0) as 
a function of the atomic resonant frequency ft for the 
coupling strength g = 0.6w by means of the numerical 
simulation, the GVM and the GRWA. In the insert of 



0.4 



g=0.6co 




0.0 0.5 1.0 1.5 

Q (in unit of w) 



FIG. 4: (Color online) The mean photon number (eJa) as a 
function of the atomic resonant frequency Q for the coupling 
strength g = 0.6u> by means of the numerical simulation (NS), 
the GVM, and the GRWA. Insert: The mean photon number 
(a' a} as a function of the coupling strength g for the atomic 
resonant frequency Q = 1.5a;. 



Fig. 4, the mean photon number (a^a) as a function of 
the coupling strength g is also plotted when the atomic 
resonant frequency is chosen as ft — 1.5w. Both these fig- 



ures show that our obtained explicit expression in Eq.© 
agrees well with the numerical simulation. It implies that 
the required mean photon number can controlled well by 
manipulating an external flux bias (the atomic resonant 
frequency), the inductive coupling of a qubit (the cou- 
pling strength), and a transmission line resonator (the 
photon frequency) in current experimental setups of the 
ultrastrong coupling 13]. 

In summary, we have presented the GVM to analyti- 
cally obtain some interesting ground-state properties of 
the Jaynes-Cummings model with the ultrastrong cou- 
pling in all regions of the atomic resonant level including 
the negative detuning, the resonant case, and especially 
the positive detuning. However, for the excited state, 
this method is invalid since the second-order perturba- 
tion energy will be larger than the unperturbed ground- 
state energy. 
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